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We discuss the spin-dependence of the effective two-body interactions appropriate for three-body 
computations. The only reasonable choice seems to be the fine and hyperfine interactions known 
for atomic electrons interacting with the nucleus. One exception is the nucleon-nucleon interaction 
imposing a different type of symmetry. We use the two-neutron halo nucleus n Li as illustration. We 
demonstrate that models with the wrong spin-dependence are basically without predictive power. 
The Pauli forbidden core and valence states must be consistently treated. 
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Introduction. Quantum halos occur in several 
branches of physics although mostly discussed in 
molecules and nuclei Q, 0, 01! ■ These structures are 
described in terms of a few weakly bound composite 
clusters. The corresponding two-body interactions 
should in principle be derived from the basic interactions 
between the particles within the clusters. However, to 
be practical effective forces must be employed for at 
least two reasons. First the calculations simplify and 
second the inherent inaccuracy in the calculations from 
first principles is often much worse than in the few-body 
computations. Thus the accuracy and inter-dependence 
of computed observables can be drastically improved 
by use of phenomcnological effective interactions. This 
division between understanding the basic interactions 
and understanding the resulting many-body system is 
strikingly illustrated by nuclear structure which was 
described virtually independent from detailed knowledge 
of the nucleon-nucleon interaction |5(. 

In descriptions of nuclear halos effective two-body in- 
teractions are indeed always constructed and used with 
only superficial connection to the fundamental forces 
H E U H ES 113 This strategy is particularly well 
suited for the spatially extended halos where details of 
the potentials are unimportant. Only low-energy scatter- 
ing properties are crucial for the gross structures. How- 
ever, more quantities are both computed and measured 
with ever increasing accuracy 0, 0, 0, 0, 0, . At 
some point we are bound to exceed the validity limit of 
the models. Designing the models to maximum perfor- 
mance is very desirable. Their range of applicability can 
then be extended by careful choices of model space and 
form of the effective interactions. 

Few-body computations with spin-dependent effec- 
tive interactions have so far not been very abundant 
IE EJ H> 13 but this is likely to change in the near fu- 
ture. It is then important a correct starting point which 



unavoidably is the two-body interaction. The purpose 
of this report is to discuss the constraints on the spin- 
dependence of the effective two-body interactions. We 
shall test by application to three-body nuclear halo sys- 
tems where the effects are more indirect, to some extent 
even hidden, and only revealed by systematic computa- 
tions. 

Symmetry requirements. It was noted early in the his- 
tory of nuclear physics that a spin-dependence in the nu- 
cleon mean-field potential was necessary to reproduce the 
magic numbers |5j. The only rotationally invariant, par- 
ity conserving, spin-dependent quantity is the spin-orbit 
potential which is doing the job and therefore used in all 
realistic mean-field computations. In complete analogy 
the nucleon-core effective interaction for zero core-spin 
can then only have central and spin-orbit terms. 

On the other hand the nucleon-nucleon interaction has 
central, spin-spin, spin-orbit, quadratic spin-orbit, and 
tensor terms |20j. For an even-even nucleus they reduce 
to central and spin-orbit terms after applying the mean- 
field average as in all optical model computations [2l| . 
However, for non-zero core-spin all these spin-dependent 
forms are allowed in the construction of the nucleon-core 
effective potential. We shall in this report for simplicity 
omit quadratic spin-orbit and tensor interactions. The 
first is usually small and the effect of the second is re- 
produced by an orbital angular momentum dependent 
central force if the non-diagonal part is neglected or does 
not contribute. In few-body physics this is consistent 
with the common choice of Hilbert space with only a few 
opposite parity relative states, e.g. s and p-states. 

This leaves us with central, spin-spin and spin-orbit 
terms. The spin-dependence of the effective two-body 
interaction must then be a combination of the three pos- 
sible independent scalar quantities, i.e. 

as n ■ s c + b£- s n + c£- s c , (1) 
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where s n and s c are the spins of the two particles, i 
is the relative orbital angular momentum operator and 
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a,b,c are constants. The spin-symmetric combination, 

as n ■ s c + b£- (s n + a c ) , (2) 

has to be used for identical particles like nucleons [20| . 
However, an asymmetric system like the electron and the 
nucleus require the combinations 

a(£ + s n ) ■ s c + b£ ■ s n , (3) 

where "n" now refers to the electron and s c is the total 
angular momentum of the nucleus. The origin of the 
first term in eq.© is the magnetic hyperfine interaction 
arising from the nuclear and electron spins [22|. 
Other combinations are clearly also possible, e.g. 

as n ■ s c + bi ■ s n , (4) 

{£ + S n ) ■ (£ + S c ) = £■£ + S n ■ S c + I ■ (S n + S c ) ,(5) 

where eq. (J2J is more symmetric but equivalent to eq.© 
with the £ ■ £ term included in the central potential. 

For the nucleon-core interaction the absence of symme- 
try requirements is complicating the choice. One option 
could be to use eq.© as in where three-body sys- 

tems with two nucleons outside a core of finite spin are 
investigated. With appropriate parameter adjustments 
the computed spectra obtained in 0] and the results 
involving both the ground state structure and breakup 
reactions are indeed quite reasonable p3| . 

When the combinations in eq.© are used the total 
spin s = s c + s n of the neutron-core system is a conserved 
quantum number. When the core has a finite spin the 
two values of s = s c ± 1/2 decouple and separate through 
the spin-spin term. Each of these two spin values is then 
split up by the spin-orbit term according to the possible 
total two-body angular momenta j 2 = (£ + s) 2 . However 
s and j are not the usual mean-field quantum numbers, 
where every nucleon moves in an orbit characterized by 
the relative nucleon-core orbital I and nucleon total j n = 
£±1/2 angular momentum. Instead j n couples to s c to 
give the total two-body angular momentum j. 

One serious problem with eq. © arises because the nu- 
cleon angular momentum j n is not conserved like the to- 
tal spin. Therefore the usual mean-field spin-orbit part- 
ners with jn = £ ± 1/2 are inevitably mixed in both two 
and three-body systems. The essence of the difficulty is 
that the motion of the valence nucleon(s) outside the core 
is inconsistent with the (approximate mean-field) motion 
of the identical nucleons within the core. This is obvi- 
ously a problem for a few-body description in terms of 
components with given j n as needed when we want to 
bridge the gap to mean- field computations. 

The problem is much worse if one and only one of the 
mixed spin-orbit partners is Pauli forbidden by core oc- 
cupation as the orbit in 10 Li and 11 Li. This may 
be disguised by parameter adjustments, but sufficiently 
many independent observables computed from the same 
parameter set would undoubtedly reveal the problem. 
Computations of relatively few observables may be quite 
reasonable although uncontrolled and thus suspicious. 



We need to restore the mean-field single-particle to- 
tal angular momentum as a conserved quantum number. 
Therefore in subsequent calculations (see for example 
|24j) the spin-orbit interaction was changed and the spin- 
dependence were instead given by eq.©. Now the domi- 
nating spin-orbit term conserves I and j n as indicated by 
the usual notation ij n . The spin-spin term should then be 
a perturbation precisely designed to split the states with 
different total angular momenta arising from couplings 
of core (s c ) and nucleon (j n ) spins. This is a reasonable 
assumption when the single particle spin-orbit splitting 
is much larger and essentially is maintained. 

To achieve full consistency with the mean-field descrip- 
tion we have to replace eq.© with eq.©, which con- 
serves jn ■ It is interesting that this choice is appropriate 
for nucleon-core as well as for for electron-nucleus inter- 
actions. In both cases the fine- and hyperfine structure 
can apparently then be properly described as originat- 
ing from the spin-spin and spin-orbit terms, respectively. 
The difference between eqs. (JSl and © is marginal for 
the present cases of interest, since the strength of the 
spin-spin term in both cases must be rather small. How- 
ever, eq.© is clearly preferable from a conceptual point 
of view where the connection to mean-field computations 
is necessary or at least very desirable. 

It is worth emphasizing that eqs.©, ©, © and © 
are identical both for s-waves and for zero core spin. 

Comparing two-body properties. The interaction must 
in general be combinations of three terms as in eq.© 
with appropriate individual radial form factors. The rel- 
ative importance of the terms in eq.© is illustrated for 
d-waves in the upper part of figH We define a = x and 
b = \/l — x 2 such that the parameter x controls the rela- 
tive weight between the spin-spin and spin-orbit parts of 
these potentials. For x — ±1 only the spin-spin term en- 
ters, the total spin is conserved, and the state with s = 
is separated from the three degenerated s = 1 states. 
When | a; | < 1 the spin-orbit part breaks this degeneracy 
leaving four different energies as for x = where only the 
spin-orbit term contributes. 

The relative importance of the terms in eq.© is il- 
lustrated in the lower part of figd where we again use 
the parameter x as a measure of the relative weights of 
spin-spin and spin-orbit potentials. For x — (only spin- 
orbit) the eigenvalues coincide for the same value of j n - 
For finite x the spin-spin interaction removes this degen- 
eracy leaving four different eigenvalues corresponding to 
the different pairs of j n and j. 

The difference between the interactions in eqs.© and 
© is illustrated in fig|21with the potential 

V(x) = s n ■ s c + B£ ■ (s n + s c ) + x(l - B)£ ■ s c , (6) 

varying from the form in eq. © for x — to that of eq. © 
for x = 1. In both cases B is the strength of the spin- 
orbit term. We have chosen a large value of B to enhance 
the effect of the spin-orbit term. For x = the total spin 
s and total angular momentum j are conserved quantum 
numbers collectively denoted by Sj. When x increases 
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FIG. 1: Eigenvalues for d-waves and s c = s n — 1/2 of the 
potentials in eqs.Q (upper part) and 10 (lower part) as func- 
tions of x where a = x and b — vl — x 2 . We use the notation 
£^K The orbital (I), nucleon (j„) and total (j) angular mo- 
menta are defined in the text. 



s and j are not conserved quantum numbers anymore, 
and the eigenvalues evolve as shown in fig|21 For x = 1 
when we have the form in eq.© the conserved quantum 
numbers {£,j n ,j) are collected in the notation lj\ 

Let us first focus on the upper part of the figure. For 
x = 1 we see clearly the spin-orbit splitting of the d§/2 
and d 3 / 2 -states. The c£ 5 / 2 -states are lower due to the 
negative sign of the strength B. This illustrates that the 
choice in eq.(j3J) (x = 1) permits construction of low- lying 
^5/2-states without contribution from the ^3/2- This is 
not possible with eq.J2Jl (x — 0), where two of the eigen- 
values mix the d 5 / 2 and <i 3 / 2 -states, while the other two 
are pure ci 3 / 2 and pure (£5/2, respectively. Thus it is not 
possible with the choice in eq.(0) to have two low-lying 
d 5 / 2 -states without admixtures from the d 3 / 2 -state. 

Similarly, the dashed lines in fig|3 are the eigenvalues 
of the potential 



V(x) = s n ■ s c + Bt ■ s n + xl ■ s c 



(7) 



that for x = and 1 reduce to eqs.|@J and respec- 
tively. We see that although j n is not a conserved quan- 
tum number even for x = 0, the eigenvalues remain al- 
most constant all the way up to x = 1. This result is 
a consequence of the small absolute size of the spin-spin 
interaction compared to the spin-orbit term. Under this 
realistic condition Q is a good approximation to J3J|. 

In the lower part of figEl we show the same kind of 
results for p-waves and s c =3/2. We start again with the 
expression eq.© but now with B = 20. The positive 
value of B reverses the normal order of the spin-orbit 
splitting. Then for x—1, where j n and j are good quan- 
tum numbers, the eigenvalues with smallest value of j n 
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FIG. 2: Eigenvalues of the potentials eqs.jfjj (solid lines) 
and @ (dashed lines). The upper part is for d- waves 
and s c —s n =l/2, and the lower part is for p-waves and 
s c =3/2,s„=l/2. These are the waves relevant for 17 Ne and 
11 Li, respectively. When conserved, the total spin s and to- 
tal j angular momentum are collected in the label Sj and 
analogously £^ labels state with conserved orbital, particle 
"n" and total j angular momentum. B is the strength of the 
spin-orbit interaction. 



are the lowest. This is a useful tool to exclude occupation 
of the Pauli forbidden p 3 / 2 -states achieved by shifting 
them to high energies. Again we see that eq.Q (x—1) 
allows close lying, still spin-split, states with equal value 
of j n . In contrast eq.@ (x = 0) corresponds to two pure 
Jn=3/2-states and four other states which mix the P1/2 
and the p 3 / 2 orbits according to 



(8a) 
(8b) 
(8c) 
(8d) 



It is then clear that when eq.© is used, the two low- 
lying p-states necessarily contain part of the p 3 / 2 -states. 
They are forbidden by the Pauli principle for both 10 Li 
and 11 Li. The dashed lines in the lower part of fig|21show 
eigenvalues obtained from the potential eq.J7J). Again 
we observe that eq.lQJ (x = 0) is a good approximation 
to eq.@ (x = 1) provided the spin-spin term only is a 
perturbation to the spin-orbit term. 

Comparing three-body properties. Let us now inves- 
tigate how the different choices of the nucleon-core in- 
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teraction affect the full three-body calculation. We take 
11 Li ( 9 Li+n+n) as an example with the Hilbert space 
consisting of s and p- waves. Then s c =3/2 and the low- 
est neutron Si/ 2 and p 3/ / 2 -orbits are fully occupied in the 
9 Li core, and therefore not available for the valence neu- 
trons due to the Pauli principle. The calculations are 
performed by solving the Faddeev equations with the 
hyperspheric adiabatic expansion method p5|. For the 
neutron-neutron interaction we use eq. (J2J) supplemented 
by a central part. The radial form factors are gaussians, 
and their parameters are adjusted to reproduce low en- 
ergy nucleon-nucleon scattering data. 

For the nucleon-core interaction we first use eq. plus 
a central part, i.e. 
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We use gaussian radial form factors adjusted indepen- 
dently for s and p-waves. To account for the Pauli prin- 
ciple for the S1/2 and p 3 / 2 states we start by using a 
simple shallow s-wave interaction without bound states 
and a large repulsive spin-orbit strength shifting the P3/2- 
states to high energy as seen in the lower part of fig for 

35=1. 

As in [2^| we choose a range for the 9 Li-neutron in- 
teraction of 2.55 fm, and adjust the strengths of Vc* 
and Vss to place the two low-lying spin-split P1/2- 
resonances in 10 Li at 0.3 and 0.5 MeV computed as 
poles of the S'-matrix. These values are consistent with 
the available experimental data [23. We finally use the 
strength of the central s-wave interaction to fit the exper- 
imental 11 Li separation energy, and the strength of the 
s-wave spin-spin potential to place a low-lying virtual s- 
state in 10 Li at 50 keV, as indicated by experiments. 

The precise values of the parameters in the potentials 
are given in the caption to figGO The computed r.m.s. 
radius is 3.3 fm, and the p-wave content is slightly above 
30%. Use of a weaker central s-wave interaction can in- 
crease the p-wave content up to around 40% in somewhat 
better agreement with the known experimental values. In 
this case an additional effective three-body potential is 
needed to maintain the two-neutron separation energy at 
the right value. However, this fine tuning is not necessary 
for the present investigation. The 10 Li energy distribu- 
tion is approximately obtained as the Fourier transform 
^(Px,Py) of the n Li wave function, i.e. 
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where p x is the 9 Li-neutron relative momentum whose 
direction is given by r2 Px , fj, x is the 9 Li-neutron reduced 
mass, p y is the relative momentum between the second 
neutron and the 10 Li center of mass, and E x = p^/2fj, x . 

Let us carry out the same calculations but now using 
eq.®. The 9 Li-neutron interaction is given as 

Vg = V W + Vjps n ■ s c + V&t ■ (s n + s c ). (11) 
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FIG. 3: The dominating radial component for the ground 
state of 11 Li (outer part) and the 10 Li energy distribution 
as given in eq. (HOP (inner part). The potentials in eqs.@ 
(solid line) and 1111 (dashed line) are used. The strengths of 



Ai=o) T/ (fco) T/ (fci) 



-7.05 MeV, -1.6 MeV, 260.25 MeV, 1.0 MeV and 300 MeV, 
respectively, for the potential eq.© and —5.95 MeV, —1.6 
MeV, -27.125 MeV, 1.5 MeV and 5 MeV for the potential 
eg, lllll . The range of the gaussians is 2.55 fm in all the cases. 



The p-resonances are now characterized by the quantum 
numbers given in the lower part of fig|21 for x=0. We 
proceed as before, placing the two lowest p-resonances 
at the energies of 0.3 MeV and 0.5 MeV and using the 
spin-orbit strength to shift the remaining four states to 
higher energies. Again the s-wave interaction is used to 
adjust to the experimental separation energy of 0.3 MeV. 
This procedure results in a r.m.s. radius of 3.4 fm, also 
consistent with the experimental data, but the p-wave 
content is now too small, less than 7%. 

An increase of this p-wave content without changing 
the energies of the two lowest resonances can only be 
achieved by reducing the strength of the spin-orbit in- 
teraction. With the parameters given in the caption of 
fig|2|we then obtain a reasonable n Li wave function, with 
the right separation energy, r.m.s. radius, and a p-wave 
content of 27%. Furthermore, as shown by the dashed 
curves in figEl the dominating radial wave function as 
well as the 10 Li energy distribution obtained with the 
potential ea. (|llfl resembles those obtained with eq.JHJ- 

However, some hidden differences remain. First, the 
p 3 / 2 -occupation probability only amounts to 7% with the 
potential in eq.® and around 13%, more than half of 
the total p-wave contribution, with the potential eq. (|ll|l . 
Second, by increasing the repulsion in the spin-orbit in- 
teraction in potential eq.@ this p 3 / 2 -probability can be 
reduced without significant changes of the computed n Li 
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properties. However, with potential ea. i|ll|) it is not pos- 
sible to reduce the p 3 / 2 -probability and still keep the p- 
wave content at a realistic fairly high value. Third with 
the potential in eq.|0 the computed 10 Li has only two 
low-lying resonances of P1/2 character, while for ea. i|ll|) 
almost all the p-resonances are rather close-lying, i.e. the 
five lowest energies (computed as poles of the ^-matrix) 
are 0.3, 0.5, 1.0, 1.3 and 1.6 MeV. 

Let us now maintain the realistic properties of the n Li 
ground state and turn to the lowest 1 _ excitation, i.e. 
the l/2 + state. The only difference in the calculation is 
that the inclusion of the Faddeev components must be 
consistent with the new quantum numbers. The same 
potentials as for the ground state must be used, since 
precisely the same partial waves constitute both ground 
state and excited state configurations. 

Now drastic differences appear. For the potential in 
eq.@ the computed l/2 + state is unbound. Further- 
more, a calculation with the complex scaling method re- 
veals a l/2 + resonance at 0.4 MeV with a width of 0.3 
MeV. These values are consistent with the detailed cal- 
culations shown in 13]. However, when the potential 
in ea. (|llf) is used the l/2 + state is bound by 1.3 MeV, 
i.e. even more bound than the ground state. We can then 
conclude that when only ground state properties are used 
to fit the potential parameters, the excited states are not 
automatically also correct even when they are spanned 
by the same Hilbert space. 

For the 1~ excitations of 11 Li basically one s and one 
p-state are simultaneously occupied unlike the ground 
state where both neutrons simultaneously are located in 
pairs of either s or p-states. The l/2 + resonance ob- 
tained with eq.© does not contain contributions from 
the p 3 / 2 -waves when a realistic spin-orbit splitting is ap- 
plied. In contrast this resonance obtained with ea. i|ll|) 
contains a p 3 / 2 -contribution of about 10%. The wrong 
properties of the interactions may be parameterized away 
while still reproducing the ground state observables, but 
catastrophically wrong results may return for the other 
properties for example of l - excitations. 

The Pauli principle. The calculations shown in the 
previous section contain an important simplification con- 
cerning the way the Pauli principle is taken into account. 
This has been done to make the conclusions clear and 
avoid the mixing with technical details that could easily 
obscure them. We now shall investigate if a better treat- 
ment of the Pauli principle, necessary for any realistic 
calculation, is modifying some of the previous results. In 
|27| an appropriate method to treat the Pauli principle 
in few-body calculations is described. The realistic two- 
body interactions able to bind a nucleon into a Pauli for- 
bidden state are substituted by the corresponding phase 
equivalent potentials with exactly the same phase shifts 
but without the Pauli forbidden bound states. 

In (2(| constraints on the neutron- 9 Li were derived con- 
sistent with the known 10 Li spectrum and the n Li prop- 
erties. In particular it was found that the interaction in 
eq. @ with gaussian radial form factors with range equal 



to 2 fin, and strengths for V^°\ V s { t°\ V c (fcl) , V# =1) , 
and V s ( o =1) equal to -94.0 MeV, -11.4 MeV, -79.64 
MeV, 1.1 MeV and -13.12 MeV, respectively, results in 
a 10 Li spectrum with a low- lying 2~ S\/2 virtual state at 
50 keV, and a 1+/2+ p 1/2 -doublet at 0.25 MeV and 0.54 
MeV, respectively, consistent with the available experi- 
mental data. This neutron- 9 Li interaction has a deeply 
bound si/2 state and a bound p 3 / 2 state at —4.1 MeV. 
These two states are Pauli forbidden and subsequently 
removed in the three-body calculation by use of the cor- 
responding phase equivalent potentials [26j. In this way 
the S\/2 and the p 3 / 2 forbidden states are explicitly ex- 
cluded from the calculation. 

As shown in |26J , the use of this two-body interaction 
results in a 11 Li ground state with a two-neutron sepa- 
ration energy of 0.30 MeV, a p-wave content of around 
40%, and a r.m.s. radius of 3.2 fm. These results are 
obtained after a fine tuning with an effective three-body 
interaction that is known to be necessary in three-body 
calculations to account for the polarizations of the par- 
ticles that are beyond those described by the two-body 
interactions. At the same time the n Li wave function ob- 
tained is consistent with the experimental invariant mass 
spectrum, core momentum distribution and angular dis- 
tribution obtained after fragmentation of 11 Li. Further- 
more, in 0| we show that the same interaction gives rise 
to a l/2 + n Li resonance at an energy above threshold 
of about 0.6 MeV and a width of 0.5 MeV. These results 
are similar to the ones obtained with the calculations in 
the previous section for the interaction given in eq.@. 

Let us perform now the same kind of calculations for 
potential (jllfl . We then use the s-wave neutron- 9 Li inter- 
action specified above, with a deeply bound si/ 2 state, 
forbidden by the Pauli principle, that is removed from 
the three-body calculation by use of the phase equivalent 
potential. For the p-wave interaction the good quantum 
numbers of the spin operators in the two-body potential 
are the total two-body spin s and the total two-body an- 
gular momentum j, and therefore the eigenfunctions mix 
the p 1 / 2 and p 3 / 2 states (see eas. ljSal) to <\8d\i ). It is then 
impossible to find an interaction binding pure p 3 / 2 states 
into a Pauli forbidden state, and with two pure low-lying 
Pi/2 resonances. This is due to the inconsistency between 
the good quantum numbers of the interaction lllfl and 
the mean-field description used for the core neutrons. 

This inconsistency cannot be cured by an appropriate 
treatment of the p-wave states. Let us nevertheless de- 
scribe various attempts. First, we take the same shallow 
p-wave potential as in the previous section (the strengths 
of the gaussians of the interaction i|ll|) are given in the 
caption of figOJ) . Proceeding in this way we observe that 
an appropriate treatment of the Pauli principle only in 
the s-waves while keeping unchanged the p-interaction is 
not giving new results compared to the previous ones. As 
before, the n Li ground state wave function can be con- 
sidered to be reasonable, but the same interaction gives 
rise to a bound l/2 + state with 1 MeV separation energy, 
clearly more bound than the experimental ground state. 
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A second option is to construct an interaction 
such that four of the p-wave eigenfunctions are bound, 
while the other two are low-lying resonances. The bound 
states can be interpreted as the four neutrons occupy- 
ing the p-shell in the 9 Li core. These four states should 
be removed from the calculation, and the remaining two 
states would correspond to the p-resonances at 0.25 MeV 
and 0.54 MeV. As mentioned this is in itself inconsistent, 
since the excluded states necessarily contain some P1/2 
contribution, and the two low-lying p-resonances con- 
tain part of the p 3 / 2 waves. In any case this can be 
achieved by taking gaussian radial potentials in ca. lll|) 

with a range of 2 fm and strengths for V^ -1 "* , Vss , and 
V s { o =1) equal to -77.6 MeV, -2.3 MeV, and 6.8 MeV. 
When this two-body p-potential is used, and the four 
bound p-states are removed from the three-body calcu- 
lation, the ground state n Li wave function needs a very 
large three-body effective interaction to recover the ex- 
perimental two-neutron separation energy of 0.3 MeV. 
Furthermore the p-wave content is extremely large (80%), 
and the r.m.s. radius (2.4 fm) is far below the experi- 
mental value. We see that in this case even the ground 
state wave function shows important deficiencies. Also 
the l/2 + state does not show any low lying resonance. 

A final attempt can be to use a repulsive core only 
for the P3/2 waves. For this purpose we can employ the 
phase equivalent potential constructed to be used with 
interaction @. However, this repulsive potential is im- 
plicitly assuming the shell model quantum numbers for 
the neutron state, and when used in combination with 
the spin operators in ea. lfTT|) the energies of the differ- 



ent two-body p-states bear no resemblance to the initial 
p-states, and some of them are even bound. 

We have now accounted for the Pauli principle in var- 
ious ways, i.e. using a shallow potential without bound 
states, using phase equivalent potentials and repulsive 
cores. The conclusion remains that the spin-dependent 
interaction must be simultaneously consistent with the 
treatment of both core and valence neutrons. No param- 
eter variations alter this conclusion. 

Conclusions. The first step towards a reliable descrip- 
tion of few-body systems is to construct appropriate two- 
body interactions between the pairs of particles. We 
discuss the spin-dependence of this effective interaction. 
Three different terms combining spin and orbital angular 
momentum are possible with the necessary rotational in- 
variance and parity conservation. Various combinations 
correspond to applications in different fields of physics. 
We focus especially on the forms appropriate for the 
nucleon-nucleon interaction and the fine and hyperfine 
interactions for electrons in the atom. 

The properties of different two-body interactions are 
first compared. Then the effects on computed three-body 
structures are investigated for n Li (n + n+ 9 Li). We 
show that parameters adjusted to reproduce the three- 
body ground state properties can lead to wrong results 
for the excited states depending on the form of the spin- 
dependent effective interactions. We conclude that only 
one choice is consistent with the mean-field treatment 
implicitly assumed for the core nucleons. Violating this 
condition compromises the treatment of the Pauli forbid- 
den states and reliable predictions are not possible. 
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